Dynamical Heterogeneity and Nonlinear Susceptibility in Short-Ranged Attractive 

Supercooled Liquids 
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Recent work has demonstrated the strong qualitative differences between the dynamics near a 
glass transition driven by short-ranged repulsion and one governed by short-ranged attraction. Here, 
we study in detail the behavior of non-linear, higher-order correlation functions that measure the 
growth of length scales associated with dynamical heterogeneity in both types of systems. We find 
that this measure is qualitatively different in the repulsive and attractive cases with regards to 
the wave vector dependence as well as the time dependence of the standard non-linear four-point 
dynamical susceptibility. We discuss the implications of these results for the general understanding 
of dynamical heterogeneity in glass-forming liquids. 
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The underlying reasons for the dramatic increase in the 
viscosity of glass-forming liquids are not well understood. 
It has become increasingly clear that simple structural 
measures remain short-ranged close to the glass transi- 
tion, and thus a growing simple static length scale does 
not appear to be implicated This has led to the 
search for a growing dynamical length scale that drives 
vitrification. Indeed^recent simulations @, H, 0> H> @> 
and experiments [1, 0, [ll , 12, HI have given direct 
evidence for both a growing length scale and a dynam- 
ical scaling relating its growth to the rapidly increasing 
time scales that characterize the glass transition. The 
study of this key aspect of dynamical heterogeneity, as 
encoded in various multi-point dynamical susceptibilities, 
has opened up the ability both to extract absolute length 
scales associated with cooperative relaxation in glassy 
systems and to provide precise metrics for the testing of 
various theoretical approaches 0, Q, EH . 

The simplest model system that exhibits the ex- 
pected dynamical behavior associated with more com- 
plicated glassy systems is the hard-sphere liquid. Here, 
entropy-driven crowding effects give rise to a character- 
istic dynamical behavior that includes a two-step non- 
exponential relaxation, a dramatic increase in relaxation 
times associated with small changes in volume fraction, 
and dynamical heterogeneity accompanied by a growing 
dynamical length scale 14j . Recentl y, it has been demon- 
strated via theory [l|| , simulation [PA [l8| , and experi- 
ment flij , [20l | that another extreme glassy limit exists 
for simple spherical particles: that of the short-range at- 
tractive glassy state. Here, strong short-ranged bonding 
between the particles can lead to extremely slow relax- 



ation, but with dramatically different dynamical charac- 
teristics. Dynamical heterogeneity also exists close to the 
attractive glassy state [13, El H3| , but has not been 
systematically characterized, and multi-point dynamical 
susceptibilities for such systems have not been measured 
or computed. The goal of this work is to investigate in 
detail the properties of standard non-linear spatiotempo- 
ral susceptibilities at distinct points along the attractive 
glass line and to quantitatively and qualitatively com- 
pare the observed features to those of the hard-sphere 
system that lies at the infinite temperature limit of the 
dynamical arrest line. The outcome of this exercise is a 
greater understanding of the physics of attractive glass- 
forming systems and dynamical heterogeneity in glass- 
forming liquids in general. Further, this work provides 
important benchmarks for the testing of various theoret- 
ical approaches. 

The systems that we consider have potentials of the 
form 



U(r) = 4e 



in 
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where the temperature scale T is set by e, the length scale 
is set by ubb and time t is rescaled by (e/rncr 2 ^) 1 / 2 . For 
our study n — 40 and 30, yielding a potential with an 
attractive range of approximately 3% and 4% of <jbb, 
respectively. To prevent crystallization, a 50:50 mix- 
ture with size ratio <jaaI&bb — 1.2 and <jab = &ba = 
(&aa + cbb)/2 is used. Standard molecular dynamics 
simulations with a number of particles N = 256 have 
been performed in the microcanonicai ensemble with a 
time step At < 1.3 x 10~ 3 . Finite-size effects were 
tested by comparison to a N = 2048 system with lit- 
tle discernable difference found for the quantities stud- 
ied here. In Fig. [1] we plot a dynamical phase diagram 
(T vs. volume fraction (f) of the system. The arrest 
line is determined by extrapolation of the iso-diffusion 
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FIG. 1: (Color online) Dynamical phase diagram for interac- 
tion ranges n — 30 (diamonds) and n — 40 (squares) obtained 
by extrapolating to the limit of zero diffusivity [24| . The lines 
are guides for the eye. The putative higher-order (Aa) sin- 
gularity lies near the n = 30 cusp. Simulation points B and 
C are indicated, while point A lies off the plot in the high-T 
limit of the T-(j> diagram. Dashed lines are guides for the eye. 




FIG. 2: Maximal values of Xi{k,i) as a function of wave vec- 
tor for points A (diamonds) and C (circles), relative to the 
absolute maximum of Xi(k,t) for that state point. The k- 
dependence of point B is qualitatively similar to that of point 
C and is not shown. The nearest-neighbor peak of S(k) for 
both state points is contained within the shaded area. Dashed 
lines are guides for the eye. 



curves to the limit of zero diffusivity [24[ . Three (T, <fi,n) 
state points in the supercooled-liquid regime are con- 
sidered: A = (4.4,0.605,30), B = (0.36,0.59,30), and 
C = (0.34,0.6,40). Point A lies close to the hard-sphere 
limit, while point C lies close to the attraction-driven ar- 
rest line, but away from the putative (A3) dynamical sin- 
gularity predicted by mode-coupling theory (MCT) [25[ . 
Point B lies close to the arrest line in the "reentrant 
pocket", near the location of the higher-order (A4) sin- 
gularity predicted by MCT. In all cases, state points 
have been chosen not only to reflect potentially distinct 
physics, but so that the a-relaxation time T a and the bulk 
diffusion constants are similar as well . This allows for 
comparison of the potentially different physics at points 
A-C for comparable degrees of absolute sluggishness. 



Simple two-point dynamical correlation functions have 
already been extensively characterized in these sys- 
tems [24|. The two-point function F s (k,t) displays dras- 
tically different decay characteristics at the three se- 
lected state points, ranging from standard two-step re- 
laxation with associated power-law relaxation at point A 
to intermediate-time logarithmic decay spanning several 
time decades at point B. Here we characterize the fluctu- 
ations of two-point dynamical quantities that are directly 
relevant for understanding dynamical heterogeneity. In 
particular, we focus on the normalized Xi{k,t) suscepti- 
bility defined as 



X4(M) 



(f s (k,tf) - (/,(fc,t)) 2 

w _1 £i [</i(M) 2 > - (/i(M)> 5 



(2) 



where f s (q,t) = N' 1 £. cos {k • [n(t) - r,(0)]} = 
N -1 J2i fl{k, t). This quantity measures the size of fluc- 
tuations in self-density correlations at a particular wave 
vector; its growth is related to the increase of a coopera- 
tive dynamical length scale. Here, unlike in some earlier 
work, we explicitly label the wave vector dependence of 
this susceptibility, which indicates that the fluctuations 
associated with dynamical heterogeneity may be large 
or small depending on the intrinsic length scale that is 
probed. Indeed, as discussed by Chandler et al. [lj| (see 
also Ref. [i3|)i tne fc-dependence of %4(fc,i) is a useful 
way to probe the various length scales associated with 
dynamical heterogeneity. The ^-dependence of Xi(k,t) 
should not be confused with that of S^k.t), which pro- 
vides a four-point analog to the static structure factor, 
and allows for the direct extraction of a length scale as- 
sociated with cooperative heterogeneous motion. On the 
other hand, it is expected that max t [xi{k, t)] ~ £(t) 2 ~ v , 
where £(t) is the same dynamical heterogeneity length 
scale extracted from S^k^t) [28| and n is the suscepti- 
bility exponent. Thus, as long as n does not vary for 
the region of interest in the dynamical phase diagram, 
one may infer some information concerning the growth 
and absolute size of £ (t) [2§| . Lastly, it should be noted 
that the definition of X4(k,t) given above slightly differs 
from the standard definition, due to the normalization 
factor in the denominator. This normalization is used 
to attempt an unbiased comparison of peak amplitudes. 
We have checked that the conclusions drawn from the 
results presented below are not altered if the standard, 
unnormalized definition is used. 

We start with a comparison between the fc-dependence 
of the maximal values of X4(k,t) for state points A and 
C. As shown in Fig. [21 a striking qualitative distinction 
exists between the size of dynamical fluctuations in the 
cases where glassy behavior is driven by strong, short- 
ranged bonding compared to the hard-sphere limit, where 
crowding drives vitrification. In particular, dynamical 
fluctuations are maximized for wave vectors below that 
of the main diffraction peak of S(k) in the hard-sphere 
limit, while in the attraction-driven case the maximal 
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fluctuations occur for wave vectors in excess of the first- 
neighbor peak of S(k). This finding makes clear the fact 
that, while in the hard-sphere case dynamical heterogene- 
ity fluctuations are most sensitive to collective events on 
scales larger than the particle size, in the case of strong 
short-ranged attractions it is bonding fluctuations that 
trigger the emergence of dynamical heterogeneity (poten- 
tially associated with large length scales) , as measured in 
Xi(k,t). 

It is instructive to compare this result with the recent 
calculations of Greenall et al. [3(| • In this work, the sensi- 
tivity of the fc-dependent plateau height to changes in the 
structure, as computed by MCT, was measured for both 
the hard-sphere and the attractive glass-forming limits. 
For the hard-sphere system, it was found that sensitivity 
is most pronounced for changes in structure just beyond 
the first shell of neighbors. Greenall et al. deemed this 
the "caged-cage" effect. On the other hand, it was found 
that the plateau for systems near the attraction-driven 
arrest line is most sensitive to changes of structure that 
occur at high wave vectors associated with short-ranged 
bonding - i.e. at fc values much in excess of an inverse 
particle size. 

A strong qualitative similarity is thus seen between 
the fc-dependence of the plateau sensitivity, as computed 
by MCT, and the fc-dependence in the peak height of 
X4(k,t) as measured directly via MD simulation. To 
interpret this we first remark that, at least in the (3- 
relaxation regime, %4(fc,i) is an indicator of the spatial 
fluctuations of the plateau height, which has been shown 
to be spatially heterogeneous even at these short time 
scales in recent numerical simulations [3l| . Clearly the 
MCT calculations of Greenall et al. measure the sen- 
sitivity of the plateau of a uniform system to uniform 
changes in structure, while at any instant in a real liq- 
uid the local structure varies from site to site. However, 
it is reasonable to assume that these spatial fluctuations 
will mirror the very same sensitivity to local structure 
as the global plateau does to a global change in struc- 
ture. Thus, our dynamical results for the fc-dependence 
of the peak height of fluctuations associated with dy- 
namical heterogeneity provides a deep connection with 
the static MCT calculations of Greenall et al. This in- 
terpretation is in harmony with recent calculations and 
speculations concerning the nature and interpretation of 
dynamical heterogeneity within MCT 32, 331 ■ 

We now turn to the full time dependence of X4(fc, t) . In 
Fig. [3] we show x±{k, t) for k- values above and below the 
peak of <S'(fc) for the three different state points. Clearly 
the temporal shape of X4(fc, t) is qualitatively different at 
the three points. The time-dependent growth of %4 (&;,£) 
to its peak in the repulsion-driven limit of point A may 
be fit to a power-law form, as shown in the log-log inset 
of Fig. [3]^top) . This behavior is similar to that observed 
in many other systems, such as mixtures of soft-sphere 
or Lennard- Jones particles [f| . This is in contrast to the 
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FIG. 3: (Color online) Time evolution of Xi(k,t) f° r wave 
vector k = 6.6 (squares), 11.2 (diamonds), 18.7 (circles), and 
22.5 (triangles) at state points A, B, and C, from top to 
bottom. Insets: thick black lines show approximate power- 
law growth at intermediate times in hard-sphere limit (top) 
and logarithmic growth at state point B (middle). Dashed 
lines are guides for the eye. 



case of point B, which lies closest to a putative MCT 
higher-order singularity. Here, as shown in Fig. [3j middle) 
the amplitude of Xi{k,t) in the /3 regime is extremely 
small, while the peak heights in the a regime are sizable 
and actually exceed those calculated in the hard-sphere 
case. Thus, the f3 regime at this state point is local in its 
physics, displaying none of the hallmarks of cooperativity 
that have already set in at short times in the hard-sphere 
case. Further, the dramatic change in temporal behavior 
of Xi(k,t) from intermediate to long times suggest the 
possibility of different length scales governing the (3 and 
a regimes, respectively, in contrast to what is usually ob- 
served in typical glassy systems 34|, |35( . Lastly, at state 
point C, where attractions dominate relaxation but far 
from the location of the reentrant elbow of the arrest line, 
a behavior with mixed properties is seen. In particular, 
the growth of X4(k, t) does not display the drastic differ- 
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ence between the (3 and a regimes, although the growth 
is quite slow, and peak values reach only modest ampli- 
tudes. 

A close inspection of the behavior of X4,(k,t) at state 
point B indicates that it grows in the (3 regime not as 
a power law, but essentially logarithmically in time (see 
linear-log inset of Fig. [3{ middle)). While the decay of the 
two-point function FJk, t) is known to be nearly logarith- 
mic in this regime [24j . it is not at all obvious that this 
should also be true for its fluctuations. It has recently 
be argued that the susceptibility Yr(fc, t) ~ iE^hM ma y 
serve as a mimic of Xi{k,t) [lj, [34j, |35|. By consider- 
ing F„(k,i), it is clear that the growth of Xi(k,t) to its 
peak should be a power law for standard repulsive sys- 
tems. On the other hand, the same considerations are not 
informative for the attractive regime, where the leading 
term of xt(&, t) in the (3 regime is not dependent on time 
at all. Indeed, argumentation based on the susceptibility 
Xr{k, t) would suggest that X4(k, t) grows logarithmically 
in time only if subleading terms in the expansion of the 
two-point function is a power series in the logarithm, as 
suggested from MOT analysis near higher-order singular- 
ities (2||. Indeed, the result presented here for Xi(k,t) 
may be taken as indirect evidence for the reality of these 
subleading terms. 

In conclusion we have systematically studied how the 
dynamical heterogeneity indicator Xi{k,t) varies along 
the arrest line in attractive colloidal systems. The 
behavior of x±(h,t) in the attraction-dominated limit 
markedly differs from that previously observed in stan- 
dard repulsion-dominated systems. First, the scale of 
fluctuations is maximized for intrinsic length scales sig- 
nificantly smaller than a particle diameter, as opposed to 
the hard-sphere case where fluctuations are maximized 
at length scales in excess of the particle size. This result 
suggests that short-ranged bonding fluctuations trigger 
dynamical heterogeneity in attractive systems, while in- 
trinsic dynamics at scales larger than the cage scale cou- 
ple most strongly to dynamical heterogeneity in repulsive 
systems. In addition, the time dependence of X4,(k,t) 
varies dramatically from one limit to the other. These 
suggest marked differences in the degree of cooperativity 
seen in attractive and repulsive cases. In particular, the 
amplitude of Xi(k,t) is much smaller in the (3 regime in 
the attractive case, and the growth of x 4 (fc,£) is loga- 
rithmically slow near the onset of reentrance, as opposed 
to the more common power-law behavior. These results 
deepen our understanding of the physics of dynamical 
heterogeneity as well as provide testable targets for theo- 
retical approaches. In particular, it would be interesting 
to apply the recent extension of MCT of Ref. li^ , which 
has successfully predicted the behavior of x<i(k, t) and the 
growth of the dynamical length scale in standard systems 
to the case of attractive glass-forming systems. 
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